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HERA ep collider provides unique information on the proton structure. High center of mass energy s = 320 GeV 
gives access to both the low Bjorken- ir domain and regime of high momentum transfers Q. Recently the HI 
collaboration reported a high precision measurement of the structure function F2 at low x leading to tight 
constraints on the sea quark and gluon densities. Both the HI and ZEUS collaborations measure the structure 
function Fl which provides an important cross check of the conventional QCD picture. This measurement is 
recently extended by HI to low Q 2 where small x corrections may play important role. An ultimate precision 
of the deep inelastic scattering cross section measurement is achieved by combining the measurements of the 
HI and ZEUS collaborations. The combined data are used as a sole input to a QCD fit to obtained HERA 
PDF set. New measurements of inclusive e~p neutral and charged current scattering cross sections by the 
ZEUS collaboration at high Q 2 improve precision in this kinematic domain. HI analysis of the DIS high Pt 
jet production cross section is used for a determination of the strong coupling constant c?g. Separation of the 
strange quark density from the total sea is obtained by the HERMES collaboration using tagged production. 
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1. Introduction 

Deep inelastic lepton-hadron scattering (DIS) is im- 
portant for the understanding of the structure of the 
proton and of the dynamics of parton interactions. 
The discovery of Bjorken scaling [l[ and its viola- 
tion at fixed target experiments triggered the devel- 
opment of the theory of strong interactions, Quantum 
Chromodynamics (QCD). Significant progress in the 
exploration of strong interactions has been achieved 
at the electron-proton collider HERA. 

The high center-of-mass energy of the ep scattering 
at HERA leads to a wide kinematic range extending 
to large values of the absolute of the four-momentum 
transfer squared, Q 2 , and to very small values of the 
Bjorken x variable. At the HERA beam energies of 
E e = 27.6 GeV for the electron and E p = 920 GeV for 
the proton, Bjorken x values as small as 10 (10 -6 ) 
are accessible for Q 2 of 10 GeV 2 (0.1 GeV 2 ) and Q 2 
values as high as 30000 GeV 2 for high x. 

HERA operation spanned over 15 years, from 1992 
until 2007, with a shutdown in 2000 - 2002 to upgrade 
luminosity and instal spin rotators, to enable longitu- 
dinal beam polarization for the colliding experiments. 
The two colliding experiments, HI and ZEUS, col- 
lected the DIS data for the whole HERA running pe- 
riod. The data before the luminosity upgrade is used 
for high precision measurements at low Q 2 and low x, 
where the scattering cross section is high. The data 
after the luminosity upgrade is focused more on high 
Q 2 analyzes and polarization studies. For the last 
three months of its operation HERA run at reduced 
proton beam energy to measure the proton structure 
function F L . 

This paper is organized as follows. The DIS 
cross-section formulae and kinematic reconstruction 
at HERA are discussed in Sec [5J The resent measure- 
ments of the DIS cross section at low Q 2 and QCD 
analysis of these data performed by the HI collabora- 
tion are presented in Sec[3J An ultimate experimental 



precision is achieved by combination of the HI and 
ZEUS data which is presented in Sec [4] These com- 
bined data are used as input for a QCD fit, termed 
HERAPDF0.2. The measurement of the structure 
function F L performed by both HI and ZEUS collab- 
oration is presented in Sec [5] The new measurement 
of the neutral and charge current scattering cross- 
sections are given in Sec 15.11 and Sec 15.21 HI mea- 
surement of the DIS jet-production cross-section and 
determination of the strong coupling constant as is 
described in Sec [51 The flavor decomposition of the 
parton densities using data collected by HI and ZEUS 
is given in Sec [7] 



2. DIS Cross Section 

The neutral current deep inelastic e p scattering 
cross section, at tree level, is given by a linear combi- 
nation of generalized structure functions. For unpo- 
larized beams it can be expressed as 
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where the electromagnetic coupling, the photon prop- 
agator and a helicity factor are absorbed in the def- 
inition of a reduced cross section cr rNC , and Y± = 

1 ± (1 — y) 2 . The functions F2, Fl and xF 3 depend 
on the electroweak parameters as: 



F 2 
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Here v e and a e are the vector and axial-vector weak 
couplings of the electron to the Z boson and k z (Q 2 ) = 
Q 2 /[(Q 2 + M 2 )(4sin 2 ecos 2 e)]. At low Q 2 , the con- 
tribution of Z exchange is negligible and oy.nc = 
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Figure 1: A view of a high y event reconstructed in the HI detector. The positron and proton beam directions are 
indicated by the arrows. For the coordinate system used at HERA the z axis points in the direction of the proton beam. 
The interaction vertex is reconstructed using the hadronic final state (thin lines) and the scattered positron (thick line) 
tracks in the central tracker. The central tracker consists of (from the beam line outwards) the silicon tracker, the drift 
chambers CJC1 and CJC2, it is surrounded by the liquid argon (LAr) calorimeter. The detector operates in a solenoidal 
magnetic field of 1.16 T. The scattered positron trajectory is reconstructed in the backward silicon tracker BST and 
the CJC1. The charge of the particle is determined using the track curvature. The positron energy is measured in the 
electromagnetic part of the SpaCal calorimeter. 



F-2 — y 2 Fr,/Y + . The contribution of the term con- 
taining the structure function F^ is only significant 
for large values of y. 

In the Quark Parton Model (QPM) the structure 
function Fl is zero |3j and the other functions in 
equation^ are given as 

(F 2 ,F2 Z ,F£) = [(el,2e u v u ,vl + al)(xU + xU) 
+ (e 2 d ,2e d v d ,v 2 d + a 2 d ){xD + xD)}, 
(xF^ z ,xF z ) = 2x[(e u a u ,v u a u )(xU - xU) 

+ {e. d ad,v d ad)(xD - xD)], (3) 

where e u , e d denote the electric charge of up- or down- 
type quarks while v u<d and a u ^ d are the vector and 
axial-vector weak couplings of the up- or down-type 
quarks to the Z boson. Here xU, xD, xU and xD 
denote the sums of up-type, of down-type and of their 
anti-quark distributions, respectively. Below the b 
quark mass threshold, these sums are related to the 
quark distributions as follows 



xU — xu + xc, 
xD = xd + xs, 



xU = xu - 
xD = xd - 



xc, 
- xs . 



(4) 



charged current e^p scattering is 
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At leading order, the e + p and e charged current scat- 
tering cross-sections are 



<cc = xU +{\-y) 2 xD, 
<cc = xU+(l- yfxD. 



(7) 



Therefore the NC and CC measurements may be 
used to determine the combined sea quark distribu- 
tion functions, U and D, and the valence quark dis- 
tributions, u v and d v . A QCD analysis in the DGLAP 
formalism Q also allows the gluon momentum distri- 
bution xg in the proton to be determined from the 
scaling violations of the data. 

For NC scattering, event kinematics can be recon- 
structed using the scattered electron as well as the 
hadronic final state particles. A typical low Q 2 , high 
y event measured by the HI detector is shown in 
figure [TJ For CC scattering, kinematics is recon- 
structed using the hadronic final state. 



where xs and xc are the strange and charm quark dis- 
tributions. Assuming symmetry between sea quarks 
and anti-quarks, the valence quark distributions result 
from 

xu v = xU — xU , xd v = xD — xD. (5) 

A reduced cross section for the inclusive unpolarized 



3. Measurements of DIS Cross Section at 
low Q 2 by H1 

Recently the HI collaboration reported new mea- 
surements of the NC cross section at low 0.2 < Q 2 < 
150 GeV 2 [!, based on data collected in 1999- 
2000. This measurement comprises several analysis. 
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Figure 2: Measurement of the structure function F2 at 
fixed Q 2 as a function of x by the HI collaboration.. The 
error bars represent the total measurement uncertainties. 
The curve represents the HI PDF2009 QCD fit. 
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Figure 3: Measurement of the structure function F2 as 
a function of Q 2 at various values of x by the HI collab- 
oration. The error bars represent the total measurement 
uncertainties. The solid curve represents the H1PDF 2009 
QCD fit for Q 2 > 3.5 GeV 2 , which is also shown extrapo- 
lated down to Q 2 = 1.5 GeV 2 (dashed). 



The lowest Q 2 < 1.5 GeV 2 is reached using a dedi- 
cated HERA run in which the interaction vertex was 
shifted in the forward direction to allow detection of 
the electrons scattered at large angles. For 1.5 < 
Q 2 < 12 GeV 2 the measurement is performed using 
a dedicated HI run with modified trigger conditions. 
Both measurements use the BST (HERA-I configura- 
tion) to reconstruct the scattered electron and mea- 
sure the event vertex, and employ events with initial 
state QED radiation from the electron to reach lowest 
Q 2 values. The higher range in 12 < Q 2 < 150 GeV 2 




Figure 4: Parton distributions as determined by the 
H1PDF 2009 QCD fit at Q 2 = 1.9 GeV 2 a) and at 
Q 2 = 10 GeV 2 b). In a) and c) the gluon and sea-quark 
densities are divided by a factor 20. The inner error bands 
show the experimental uncertainty, the middle error bands 
include the theoretical model uncertainties of the fit as- 
sumptions, and the outer error band represents the total 
uncertainty including the parameterization uncertainty. . 

is measured under nominal beam and detector condi- 
tions. The event vertex is reconstructed by the CJC 
and the electron angle is measured using the vertex 
and the track element reconstructed in the backward 
drift chamber, which was mounted in front of the 
SpaCal for the HERA-I run period. 

The measured structure function F2 is shown in 
figure as a function of x for 12 < Q 2 < 150 GeV 2 . 
For each Q 2 bin, there is a strong rise of the structure 
functions for x — > 0, which indicates a large sea quark 
density. The rate of the rise increases with increas- 
ing Q 2 , which is a signature of a large gluon density. 
Figure [3] shows the structure function F2 for fixed val- 
ues of x as a function of Q 2 . For low x, there is a 
strong Q 2 dependence of the structure function. For 
highest x, the data show approximate scaling behav- 
ior. In both figures [2] and [3] the data is compared to 
an next-to- leading order (NLO) QCD fit, which de- 
scribes data very well. This fit is discussed in the 
next section. 

3.1. H1 PDF 2009 QCD Fit 

The new low Q 2 HI data together with published 
previously HI measurements of NC and CC e ± p 
scattering cross sections are used as a sole in- 
put in the NLO QCD analysis, termed H1PDF2009. 
The QCD evolution is performed the QCDNUM pro- 
gram 8]. The fit uses heavy-quark variable-flavor- 
number scheme of [j| with recent modifications [Toj |. 

The parton densities are parameterized at the start- 
ing scale Ql = 1.9 GeV 2 using xf(x) = Ax B (l + 
x) c (1 + P\x + P 2 x 2 + ...) form. Five PDFs are used to 
describe the HI data: xg, xU, xD, xu v and xd v . The 
QCD analysis imposes momentum and quark-number 
sum rules, which determine A g , A Uv and . Starting 
from all polynomial coefficients P; L = 0, the central fit 
is chosen by adding one Pi for one PDF at a time. The 
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Figure 5: Measurement of the e + p NC scattering cross 
section based on the combination of the HI and ZEUS 
data compared to the fit to these data and measurements 
from fixed-target experiments. 

modified PDF is selected if (i) x 2 is significantly im- 
proved, (ii) structure functions F2 and Fl are positive, 
(iii) all PDFs are positive, (iv) PDFs satisfy valence- 
dominance condition: xu v > xu and xd v > xd at high 
x. This analysis is repeated for all i < 3. 

The QCD analysis distinguishes three sources of 
PDF uncertainties: experimental, model and parame- 
terization uncertainty. The experimental uncertainty 
is evaluated using A% 2 = 1 criterion. The model 
uncertainty is measured by modifying the input pa- 
rameters such as masses of the heavy quarks and the 
starting scale. The parameterization uncertainty is 
estimated by adding more Pj terms (i < 3) and lifting 
the conditions (iii) and (iv). The total uncertainty 
is given by a sum in quadrature of the experimental, 
model and parameterization uncertainties. 

The results of the fit are shown in figure. 2] for the 
starting scale Q 2 = 1.9 GeV 2 and for Q 2 = 10 GeV 2 . 
The experimental uncertainty is small for all PDFs 
at all values of x compared to the total uncertainty. 
At large x, the valence quark densities have signifi- 
cant parameterization uncertainty. At low x, model 
uncertainty plays bigger role. 



4. Combination of the H1 and ZEUS data 

The highest precision of the cross-sections measure- 
ments at HERA is obtained by combining the results 
from the HI and ZEUS experiments. The combination 
of the results is performed taking into account the cor- 



related systematic uncertainties |5j, [llj . The starting 
point is the % 2 function for individual measurement 
which is defined as 



Xcxp {m,b) = 

( 8 ) 

Here // is the measured value at a point i and 7] , <5i )S tat 
and <5i iUncor are relative correlated systematic, relative 
statistical and relative uncorrelated systematic uncer- 
tainties, respectively. The function Xexp depends on 
the predictions m l for the measurements (denoted as 
the vector m) and the shifts of correlated systematic 
error sources bj (denoted as b). For the reduced cross- 
section measurements \i % = tr*, i denotes a (x, Q 2 ) 
point, and the summation over j extends over all cor- 
related systematic sources. The predictions m' are 
given by the assumption that there is a single true 
value of the cross section corresponding to each data 
point i and each process, neutral or charged current 
e + p or e~p scattering. Under the assumption that the 
statistical uncertainties are proportional to the square 
root of the number of events and that the systematic 
uncertainties are proportional to m, the minimum of 
equation [5] provides an unbiased estimator of m. 

Several data sets providing a number of measure- 
ments are represented by a total \ 2 function, which is 
built from the sum of the % 2 xp functions for each data 
set e The data averaging procedure allows the rear- 
rangement of the total \ 2 such that it takes a form 
similar to equation [8] 

The averaging procedure is applied to the HI and 
ZEUS inclusive data from the HERA-I running pe- 
riod. All the NC and CC cross-section data from HI 
and ZEUS are combined in one simultaneous mini- 
mization. Therefore resulting shifts of the correlated 
systematic uncertainties propagate coherently to both 
CC and NC data. In total 1402 data points are 
combined to 741 cross-section measurements. The 
data show good consistency, with \ 2 1 ™dof = 637/656. 
In figure [5] the NC reduced cross section, for Q 2 > 
1 GeV 2 , is shown as a function of Q 2 for the HERA 
combined e + p data and for fixed-target data [HI, [l3| 
across the whole of the measured kinematic plane. 
The data precision reaches 1.1% for 10 < Q 2 < 
100 GeV 2 . 

The combined data uses as a sole input for the NLO 
QCD analysis. The fit uses similar approach as the 
H1PDF2009 fit. The result of the fit is compared to 
the data in figure [5] There is a good agreement be- 
tween the extrapolation of the fit to lower Q 2 and 
fixed-target data in this kinematic domain. 
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Figure 6: Measurement of the structure functions Fl and 
F2 by the ZEUS collaboration compared to the prediction 
of ZEUS- JETS fit. 
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Figure 7: Preliminary measurement of the structure func- 
tions Fl and F2 by the HI collaboration. The data are 
quoted at fixed values of Q 2 and x indicated in grey and 
compared to predictions based on various models. 



5. Measurements of the Structure 
Function F L 

The large gluon density, determined from scaling 
violation of the F2 data using QCD fits (see figure H|), 
implies that the structure function Fl must be signif- 



icant at low x. Direct measurements of Fl allow to 
test this prediction providing a check of the QCD and 
adding an extra constraint for the gluon density. 

To determine the two structure functions F 2 (x, Q 2 ) 
and F L (x, Q 2 ) from the reduced cross section it is nec- 
essary to perform measurements at the same values of 
x and Q 2 but different y. This is achieved at HERA 
by reducing the proton beam energy. Two e + p runs 
at reduced proton beam energy E p — 460 GeV and 
E p = 575 GeV were performed with an integrated lu- 
minosity of about 13 pb -1 and 6 pb -1 , respectively. 
The run at E p = 460 GeV gives highest sensitivity 
to Fl while the run at E p = 575 GeV extends the 
kinematic range of the measurement and provides an 
important cross check. 

The measurement must extend to as high y as pos- 
sible to increase sensitivity to Fl. A high y kine- 
matic domain at low Q 2 corresponds to low energies 
of the scattered positron E' e . Measurement at low 
E' e is challenging primarily because of high hadronic 
background. 

The ZEUS collaboration uses Monte Carlo (MC) 
simulation to estimate the background. The MC pre- 
diction is normalized to the data using a sub-sample 
of tagged background events. The results of the ZEUS 
analysis [HI are shown in figure. [6] ZEUS report mea- 
surement of both structure functions F 2 and Fl ob- 
tained from a linear fit of the reduced cross section 
as a function of y 2 /Y + . The structure function Fl is 
found to be in a good agreement with the prediction 
of the ZEUS- JETS PDF set [H[. 

To reduce and estimate the hadronic background, 
HI demand the scattered positron candidate to have 
a reconstructed track with well measured curvature. 
The curvature is used to determine the candidate 
charge. For the signal, positive charge is expected. 
The background is approximately charge symmetric, 
i.e. the number of background events with different 
charges is about equal. Using this, the background is 
estimated from the negative charge sample, corrected 
for a small charge asymmetry and then subtracted 
from the positive charge sample. The charge asymme- 
try of the background is determined directly from the 
data by comparing negatively and positively charged 
candidates from e + p and e~p data taking periods, re- 
spectively. Therefore, the background determination 
is purely data driven for the HI analysis. 

HI published the first measurement of Fl at HERA 
using drift chambers CJC1 and CJC2 together with 
the SpaCal calorimeter for 12 < Q 2 < 90 GeV 2 [HI 
and also reported a preliminary result at higher Q 2 
using the LAr calorimeter [It}- These measurements 
were recently extended to lower 2.5 < Q 2 < 12 GeV 2 
by using the backward silicon tracker BST [l8l |. 

The HI measurement of the structure function Fl 
is shown in figure [7J The data are compared to 
the predictions from various models. The predic- 
tions agree among each other and with the data for 
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Figure 8: NC e~p scattering cross section measured by 
the ZEUS collaboration. Data for positive (red stars) and 
negative (blue dots) polarization are compared to the fit. 

Q 2 > 10 GeV 2 . For lower Q 2 , there is a notable dif- 
ference between the NLO predictions of MSTW and 
CTEQ. This difference is traced down to the difference 
in accounting for a 2 s corrections. The data are some- 
what higher than both predictions and agree better 
with the CTEQ calculations. 

5.1. Neutral Current e~p DIS Cross 
Section 

Measurements at high Q 2 using data collected in 
years 2003-2007 are published by the ZEUS collab- 
oration. For the NC scattering, taking into account 
the longitudinal beam polarization and ignoring small 
contribution from pure Z exchange, the structure 
function F 2 is 

F 2 =F 2 ~ (v e - P e a e )K Z Ff , (9) 
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Figure 9: CC e~p scattering cross section measured by 
the ZEUS collaboration. Data for positive (red stars) and 
negative (blue dots) polarization are compared to the fit. 

5.2. Charged Current e ± p DIS Cross 
Section 

In the standard model, the CC cross section de- 
pends linearly on the longitudinal electron beam po- 
larization. Figures O and [TU] show published (2(| and 
preliminary analyses of e~p and e + p CC cross sections 
performed by the ZEUS collaboration. For the pub- 
lished e~p sample, the same 2005-2006 data are used 
as for the published NC sample, the preliminary e + p 
analysis is based on data collected in 2003-2004 and 
2006-2007. The data agree well with the expectations 
of the QCD fit. The e + p CC data are of additional 
value for the QCD analyses because they are sensitive 
to the d and s quark densities which are less constraint 
by the NC data, see equations [3] and [7] 



where P e is the polarization value. Therefore, there 
is a weak polarization dependence of the cross section 
which increases at high Q 2 . 

The ZEUS measurement of the NC e~p cross sec- 
tion for P e = +29% and P e = -27% QJ] is shown in 
figure [5J The result is based on integrated luminosity 
of 169.9 pb" 1 collected in 2005 and 2006. The data 
are compared to a preliminary QCD fit. At high Q 2 , 
there is an increasing difference between cross sections 
for the two polarizations. 



6. Jet Cross Section and Determination 

of a s 

The production of jets in DIS can be used to de- 
termine the gluon density and to measure the strong 
coupling constant as- Recently HI have perfomed a 
measurement of inclusive, 2-jet and 3-jet cross sec- 
tions 21] as a function of Q 2 using HERA data col- 



lected in 1999-2007 with an integrated luminosity of 
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Figure 11: The value of as(Q) obtained by HI by a simul- 
taneous fit of all normalized jet cross sections. 
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Figure 10: Preliminary measurement of CC e + p scattering 
cross by the ZEUS collaboration. Data for positive (red 
stars) and negative (blue dots) polarization are compared 
to the fit. 



395 pb^ 1 . The measurements are well described by 
NLO QCD calculations, corrected for hadronization 
effects. 

Using these data, HI determine the strong coupling 
constant as- The measurement is performed sepa- 
rately for different Q 2 bins and different processes. 
These measurements are combined together and the 
evolution of as as a function of Q is compared to 
the theory prediction in figure [TT1 A good agreement 
is observed between the data and theory. The theo- 
retical uncertainty is dominated by the higher order 
corrections. The value of as at Mz 



a s (Mz) = 0.1176 ±0.0020(exp.)+° °^6(th.) 
± 0.0016(PDF) 



(10) 



agrees well with the world average. The experimental 
error on as is about 0.6%. The total uncertainty is 
dominated by the theory, it may be improved with 
calculation of NNLO corrections. 



7. Measurements of Strange, Charm and 
Bottom Quark Densities 

Inclusive NC cross section at low Q 2 and low x are 
dominated by the structure function F2 and do not al- 



low to separate contributions of individual quark fla- 
vors. The flavor separation can be achieved in semi- 
inclusive scattering by tagging the struck quark fla- 
vor. For HERA kinematics, tagging production of a 
charmed meson, e.g. D* , almost certainly corresponds 
to a scattering off a c-quark. Samples with secondary 
vertices are enriched with c and 6-quark scatterings. 

To measure b and c structure functions using sec- 
ondary vertices, it is essential to have high preci- 
sion tracking detector installed close to the interac- 
tion point. Recently the HI collaboration reported 
the measurement of c and b reduced cross sections (22| , 
using secondary vertex technique, based on complete 
sample for which the central silicon tracker was in- 
stalled. These data are shown in Figs [T2l and [T3l and 
compared to the H1PDF2009 fit. For all Q 2 bins there 
is a strong rise of the reduced cross section to low x 
values which increases with increasing Q 2 . This is a 
direct indication of the large gluon density. The data 
are in a good agreement with the H1PDF2009 fit. The 
uncertainty of the fit is dominated by the model uncer- 
tainty due to the variation of the heavy quark masses. 
Therefore, these data allow to check the model and 
determine these parameters. 

The fixed-target experiment HERMES, which oper- 
ated using HERA beams, recently reported a mea- 
surement the momentum distribution of the strange 
quark sea, xs(x), in scattering off the deuteron. The 
strange density is extracted from K multiplicities by 
correcting for the it, d and s fragmentation. The result 
is shown in figure [14] and compared to the CTEQ6L 
expectation. The HERMES data shows softer behav- 
ior of xs(x) compared to the CTEQ6L fit: at high x 
the data is below the fit and at low x the data are 
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Figure 12: The reduced cross section a cc measured by HI 
and shown as a function of x for different Q 2 bins. The 
predictions of H1PDF 2009 and MSTW08 NLO fits to in- 
clusive data are also shown. 



above the fit. At low x the strange sea is comparable 
to the light sea density. 
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Figure 13: The reduced cross section a measured by HI 
and shown as a function of x for different Q 2 bins. The 
predictions of H1PDF 2009 and MSTW08 NLO fits to in- 
clusive data are also shown. 



8. Summary 

Recent results from HERA provide new precise data 
for determination of the proton structure at low x. 
The combination of the HI and ZEUS results reaches 
1.1% for the 10 < Q 2 < 100 GeV 2 kinematic domain. 
The data are well described by the NLO QCD fit. This 
fit has impressively small experimental uncertainties, 
the total uncertainties are dominated by the model 
and parameterization errors. 

The conventional QCD picture is checked by the 
measurements of the structure function Fl performed 
by HI and ZEUS. Good agreement between the theory 
and the measurements is observed for Q 2 > 10 GeV 2 . 
For lower 2.5 < Q 2 < 10 GeV 2 , the preliminary data 
from HI are somewhat above the expectations. 

Measurements at high Q 2 using the polarized 
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Figure 14: The strange parton distribution xs(x) from the 
measured HERMES multiplicity for charged kaons evolved 
to Q 2 = 2.5 GeV 2 . The solid curve is the HERMES fit to 
the data, the dashed curve gives xs(x) from the CTEQ6L 
set, and dot-dashed curve us the sum of light antiquarks 
from the CTEQ6L set. 
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HERA-II data performed by ZEUS check the Stan- 
dard Model and provide constraints for the parton 
densities at high x. The flavor decomposition at high 
x is achieved by using the charged current data. 

The flavor decomposition at low x and low Q 2 
is performed by measuring semi-inclusive processes: 
tagged heavy flavor production at HI and pro- 
duction at HERMES. The HI data for the charm and 
bottom quark parton densities agree well with the 
predictions from the QCD fit to the inclusive data. 
The HERMES data are softer than prediction of the 
CTEQ6L analysis. 

Further results from HERA are expected as the data 
analysis are being finalized. In particular, combined 
HI and ZEUS measurements of the inclusive HERA- 
II cross sections and combined measurement of the 
charm and bottom production will have significant im- 
pact on PDFs for the kinematic range important for 
the future measurements at the LHC. 
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